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S : ABSTRACT 

New H I synthesis data have been obtained for six face-on galaxies with the Very 
Large Array. These data and reanalyses of three additional data sets make up a sam- 
ple of nine face-on galaxies analyzed for deviations from axisymmetry in morphology 
and dynamics. This sample represents a subsample of galaxies already analyzed for 
morphological symmetry properties in the i?-band. Four quantitative measures of dy- 
namical nonaxisymmetry are compared to one another and to the quantitative measures 
of morphological asymmetry in H I and -R-band to investigate the relationships between 
nonaxisymmetric morphology and dynamics. We find no significant relationship between 
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asymmetric morphology and most of the dynamical measures in our sample. A possi- 
ble relationship is found, however, between morphology and dynamical position angle 
differences between approaching and receding sides of the galaxy. 

Subject headings: galaxies: kinematics and dynamics — galaxies: structure — methods: 
observational 
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1. Introduction 

Despite the fact that most studies of spiral galaxy 
dynamics concentrate on understanding the proper- 
ties of axisyrnrnetric disks, evidence is accumulating 
that many galaxies lack such overall symmetry. Bald- 
win et al. (1980) were first to seriously examine the 
asymmetries in galactic disks, pointing out that lop- 
sided galaxies were a common phenomenon not local- 
ized to interacting pairs. 

More recently, frequency of morphological asym- 
metry in galactic disks has been quantitatively stud- 
ied at optical wavelengths. Based on optical mor- 
phology, approximately 30% of disk galaxies exhibit 
significant "lopsidedness" (Rix & Zaritsky 1995; Ko- 
rnreich et al. 1998, hereafter KHL). These studies are 
based on data sets containing some 30 targets each. 

Dynamical asymmetry, too, has been examined in 
large samples of single dish H I line profiles. Richter 
& Sancisi (1994), followed by the 104 galaxy sample of 
Haynes et al. (1998, hereafter HHMRvZ), examined 
the symmetry properties of single dish H I line pro- 
files and determined that as many as ~ 50% of spiral 
galaxies show departures from the expected symmet- 
ric two-horned profile. Nevertheless, asymmetries in 
line profiles are ambiguous evidence at best for dis- 
turbed dynamics, since they combine both dynamic 
and spatial information. 

Until recently, however, little work had been done 
to examine the symmetry properties of neutral hy- 
drogen in synthesis data. As a result, the connection 
between disturbed morphology and disturbed dynam- 
ics in field disk galaxies is only now beginning to be 
seriously studied. Schoenmakers et al. (1997, here- 
after SFdZ) outline a method for measuring small de- 
viations from axisymmetry of the potential of a filled 
gas disk by breaking down the observed velocity field 
into its harmonic components. Recently, Swaters et 
al. (1999, hereafter S3vA) have applied this method 
to the H I synthesis data of two galaxies, DDO 9 and 
NGC 4395, where the hallmark of dynamical asymme- 
try is found to be asymmetry in the rotation curve, 
where one side of the curve rises more steeply than 
the other. H I synthesis observations of several other 
galaxies, e.g. NGC 3631 (Knapen 1997), NGC 5474 
(Rownd et al. 1994, hereafter RDH), and NGC 7217 
(Buta et al. 1995), have also curiously revealed asym- 
metries or offsets between the optical centers of light 
and the kinematic centers of the neutral gas. While 
much evidence exists that many field spirals exhibit 



non axisymmetric morphology, and that many field 
spirals exhibit asymmetric dynamics, the question of 
whether the two phenomena are correlated, or even 
represent two pictures of a single underlying effect, 
remains open. 

Both optical morphology and gas dynamics pro- 
vide clues as to the overall structure of a galaxy. For 
instance, Zaritsky & Rix (1997) proposed that op- 
tical lopsidedness arises from tidal interactions, mi- 
nor mergers, or possibly gradual accretion. Conselice 
et al. (2000), on the other hand, find a correla- 
tion between optical asymmetry and B—V color, and 
are able to use morphological asymmetry to identify 
whether starbursts in a given galaxy are likely caused 
by interactions and mergers. Similarly, the neutral 
hydrogen dynamics of asymmetric galaxies should be 
able to distinguish between tidally deformed galaxies 
and those which are dynamically isolated. For in- 
stance, while the strongly optically-lopsided galaxy 
NGC 5474 is well-known to be under the tidal influ- 
ence of its neighbor MlOl, the disturbed morpholo- 
gies of the other relatively isolated objects in the KHL 
sample are not well-explained by standard tidal inter- 
action models, which require particular dynamics as 
well as interactions with (unobserved) nearby com- 
panions. 

Alternatively, asymmetries may arise from the ex- 
citation of unstable, one-armed spiral modes possibly 
triggered by a past interaction or minor merger (Taga 
& lye 1998a, 1998b; Lovelace et al. 1999). Librations 
of the optical galaxy about the minimum of the grav- 
itational potential might also be set in motion by a 
previous interaction. As density is dependent on ra- 
dius in a galaxy, one might expect that the natural 
modes of a galaxy also depend on radius. The result- 
ing differential oscillation might result in both lop- 
sided appearance and kinematic decoupling from the 
optical light distribution. An understanding of these 
modes, if they exist, could provide a direct measure 
of the dark matter distribution (Jog 1997; SFdZ). 

Another type of nonaxisymmctry, warping of galac- 
tic disks due to non-planar motions, has been pro- 
posed as an indicator of inclined flattened halo po- 
tentials (e.g. Dckel & Shlosman 1983; Toomre 1983), 
which are required to stabilize the warp against differ- 
ential precession, as well as observational indicators of 
massive dark halos (Tubbs & Sanders 1979). Warps 
may also be due to tidal interactions and accretion, 
and have been reported to be related to non-circular 
motions, particularly m = 1 modes (Weinberg 1998). 



3 



Although "sloshing" librations in the plane of the 
galaxy would not imply a correlation between warps 
and lopsided dynamics, "flapping" librations normal 
to the galactic plane could contribute to warping. 

In this paper, we present HI synthesis data for nine 
galaxies whose optical asymmetry properties were 
quantified by KHL. Data obtained from the archives 
at the Very Large Array (VLA)|^ have been reana- 
lyzed for the galaxies NGC 5474, NGC 5701, and 
UGC 12732, and new VLA data are presented for the 
galaxies NGC 991, NGC 1024, UGC 3685, NGC 3596, 
UGC 6420, and NGC 4688. These galaxies represent 
a sample of face-on galaxies selected on the basis of 
their optical properties. 

In we discuss the acquisition and reduction of 
the new and archived data for our sample of nine 
galaxies. In §^ we analyze the data obtained for devi- 
ations from morphological and dynamical axisymme- 
try and describe methods of determining the magni- 
tudes of such deviations. These quantitative measures 
are useful for quantifying warps in face-on galaxies 
and for non-circular motions in more inclined galax- 
ies. We then use the sample to draw correlations be- 
tween the symmetry parameters of the galaxies to one 
another and to global H I properties. Finally, in 
we present our conclusions in the context of other re- 
cent work on galaxy asymmetry. 

2. Observations and Data Reduction 
2.1. Target Selection and Observations 

The galaxy sample used in this work includes nine 
galaxies from the sample of KHL. Table |l| summarizes 
the catalogued data from the Third Reference Cata- 
logue of Bright Galaxies (de Vaucouleurs et al. 1991, 
hereafter RGB) for each of the targets: column (1) 
gives the NGC or UGC designations of the targets; 
columns (2) and (3) give the right ascension and dec- 
lination of the targets in B1950.0 coordinates; column 
(4) lists the morphological types; column (5) lists the 
major and minor blue diameters, D25 and ^25; and 
column (6) gives the face-on blue apparent magnitude 
corrected for internal and galactic extinction, Bj,. 

2.1.1. Original Data Acquisition 

Of the optical sample in KHL, the six galax- 
ies NGC 991, NGC 1042, UGC 3685, NGC 3596, 



^The VLA is a facility of the National Radio Astronomy 
Observatory. 



UGC 6420, and NGC 4688 were selected for H I 21-cm 
synthesis observations based on their apparent degree 
of asymmetry and position in the sky such that they 
could be observed at night during the time sched- 
uled for the CS array. As the photometric analysis 
conducted in KHL had not yet been completed at 
the time of the target selection, estimates of asym- 
metry were made by eye. Night-time observations 
were preferred in order to minimize solar interference. 
The KHL sample was an optically selected sample of 
disk galaxies chosen on the basis of face-on appear- 
ance to permit estimation of the optical morphological 
asymmetry unencumbered by factors such as absorp- 
tion by dust and inclination errors. The KHL selec- 
tion criteria were: small axial ratio, narrow (typically 
W^Q < 100 km s~^) velocity width obtained from 
the private database of Giovanelli & Haynes known 
as the Arecibo General Catalog (AGC), and (except 
for NGC 5474 and NGC 1042) isolation from obvious 
companions. See KHL for details of the optical target 
selection. 

H I observations of these targets were conducted in 
August 1997, December 1998, and January 1999 with 
the VLA in its CS configuration. The CS configu- 
ration is a modification of the C configuration which 
provides the resolution of the C array while main- 
taining a sufficient number of short spacings to pro- 
vide sampling of large-scale structure by positioning 
one or more antennas at D configuration locations. 
In 1997, the CS configuration consisted of the stan- 
dard C array with two antennae (from the middle of 
the east and west arms) relocated at inner D config- 
uration stations. In 1998 and 1999, one antenna was 
taken from the middle of the north arm and placed at 
an inner D configuration station to provide the req- 
uisite short baselines. 

Each galaxy was observed for a total of approxi- 
mately 5.5 hours on-source integration time in 2 AD 
correlator mode with 127 channels of 12.207 kHz sep- 
aration after Hanning smoothing. This experimental 
setup provided for a total bandwidth of 1.5625 MHz 
and, at 1.4 GHz, for 2.6 km s^^ resolution. Obser- 
vations were flux and bandpass calibrated by observ- 
ing flux calibrators for 15 minutes immediately before 
and after the target observations, and phase calibra- 
tors for 5 minutes at 30 minute intervals during the 
target observations. 
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Table 1 
Summary of Galaxy Sample 













(0) 


Galaxy 


a 


6 


Type 








(B1950) 


(B1950) 


(RC3) 


(arcmin) 




Original Observations 












NGC 991 


023302.1 


-072234 


SAB(rs)c 


2.7 X 2.4 


12.30 


NGC 1042 


023756.3 


-083850 


SAB(rs)cd 


4.7 X 3.6 


11.38 


UGC 3685 


070433.1 


+614029 


SB(rs)b 


3.3 X 2.8 


12.30 


NGC 3596 


111227.9 


+150338 


SAB(rs)c 


4.0 X 3.8 


11.67 


UGC 6429 


112225.1 


+640011 


SA(rs)c 


2.1 X 1.7 


13.63 


NGC 4688 


124513.8 


+043558 


SB(s)cd 


3.2 X 2.8 


12.47 


Archive Data 












NGC 5474 


140314.6 


+535403 


SA(s)cd Pec 


4.8 X 4.3 


11.27 


NGC 5701 


143641.5 


+053450 


(R)SB(rs)0/a 


4.3 X 4.1 


11.69 


UGC 12732 


233808.8 


+255727 


Sm 


3.0 X 2.8 


13.55 



2.1.2. VLA Archival Data 

In addition to the data obtained by us at the VLA, 
several data sets obtained by other observers were re- 
covered from the VLA data archives for reanalysis. 
These data sets include observations of the galax- 
ies NGC 5474 (RDH), NGC 5701 (observed by S. E. 
Schneider in 1986), and UGC 12732 (Schulman et al. 
1997). These targets were chosen because each was 
analyzed for optical morphological asymmetry as part 
of the KHL sample, and VLA archive H I line data 
of sufficient sensitivity (about 10^" cm~^) and spec- 
tral resolution (less than 15 km s^^) for the current 
analysis existed for them. 

The archival data were retrieved as raw UV vis- 
ibility data, and reduced by the same methods as 
were used for our original observations. The raw 
data sets for NGC 5701 and UGC 12732 were Ban- 
ning smoothed on-line by their observers; that of 
NGC 5474 was not. 

An observational summary of the nine VLA data 
sets in the total sample is presented in Table ^. 
Columns in this table represent: (1) The NGC or 
UGC designation of the target; (2) The VLA con- 
figuration in which the data were obtained; (3) The 
date of the observation; (4) The observer or paper 
in which the observation was reported; (5) The to- 
tal on-source integration time, in minutes; (6) The 



total bandwidth in MHz; (7) The channel separation 
in km s~^; (8) The total number of channels in the 
data set; (9) The synthesized clean beam size in our 
natural weighted maps, in arcseconds; and (10) The 
RMS noise a observed in the line-free channels, in 
mjy beam^^. 

The archival data sets were obtained between early 
1983 and 1993, with the earlier data sets consisting 
of a small number of channels in the D array. While 
we expected to improve marginally upon the previ- 
ous data reduction for these archival data sets, we 
were surprised to find that due to the exponential 
improvement in computer technology since that time, 
the reanalysis of these data yielded a remarkable im- 
provement over the initial analyses. By cleaning the 
data to deeper levels than was possible in the past, 
we were able to measure H I column densities up 
to almost an order of magnitude fainter than those 
previously reported for the same data sets. As a 
specific example, we were able to iteratively clean 
the UV data for NGC 5474 for 200,000 iterations, 
down to the la theoretical noise level, as compared 
to the 500 iterations reported by RDH. Our reanal- 
ysis was thus able to trace the H I to a Icr column 
density of 6 x 10^* cm^^ chan^^, compared to the 
1.2 X 10^^ cm"^ chan^^ quoted in RDH, an improve- 
ment of about half an order of magnitude with the 
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Table 2 
Observational Summary 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Galaxy 


Array 


Date Obs. 


Reference 


tint 


Wband 


Wchan 


^ chart 


Beam 


a 










(min) 


(MHz) 


(km s^-'^) 




(arcsec) 


(mJy beam~^) 


NGC 991 


CS 


Aug 97 


This Work 


381 


1.5625 


2.60 


127 


25 X 20 


0.8 


NGC 1042 


CS 


Aug 97 


This Work 


380 


1.5625 


2.60 


127 


23 X 18 


0.8 


UGC 3685 


CS 


Aug 97 


This Work 


380.5 


1.5625 


2.60 


127 


19 X 18 


0.7 


NGC 3596 


CS 


Jan 99 


This Work 


503.5 


1.5625 


2.60 


127 


19 X 18 


0.7 


UGC 6429 


CS 


Jan 99 


This Work 


377 


1.5625 


2.64 


127 


20 X 16 


0.6 


NGC 4688 


CS 


Dec 98 


This Work 


370.5 


1.5625 


2.59 


127 


21 X 18 


0.8 


NGC 5474 


D 


Jul 83 


Rownd et al. (1994) 


287.5 


3.125 


2.58 


31 


54 X 51 


2.2 


NGC 5701 


D 


Apr 87 


S. E. Schneider 


750.5 


3.125 


10.41 


31 


70 X 63 


1.3 


UGC 12732 


D 


Aug 92 


Schulman et al. (1997) 


660 


3.125 


10.36 


63 


54 X 53 


0.4 



same data set. This greater sensitivity allowed us 
to uncover the very peculiar outer dynamics of this 
galaxy not observed by RDH, structure which is of 
dynamical importance. As expected, such large gains 
in sensitivity were most pronounced from the earlier 
data sets, and of less importance for later ones. 

2.2. UV Processing and Data Analysis 

Reduction of the visibility data was conducted us- 
ing the Astronomical Image Processing System (AIPS)^ 
(Napier et al. 1983). Data were examined via the pro- 
gram TVFLG, in which poor quality UV data were 
identified and flagged. Following the flagging proce- 
dure, a linear interpolation of bandpass calibrations 
obtained from the phase calibrators for each galaxy 
was obtained and applied to the target data. 

The continuum was then subtracted from each 
data set using the program UVLIN (Cornwell et al. 
1992), which subtracts continuum emission from all 
channels based on a linear fit of the visibilities in 
a given number of line-free channels, as selected by 
manual inspection. The total number of line-free 
channels used varied according to the data set, but 
typically numbered about 30 for originally acquired 
data, and 10 for archive data. The UV data were 
then transformed to the XY plane using the program 

'AIPS is distributed by the National Radio Astronomy 
Observatory. 



IMAGR. 

At this point in the calibration, it was discovered 
that the catalogued declination of the phase calibrator 
used with the UGC 6429 data set was incorrect by one 
arcminute. This had the effect of producing incorrect 
coordinates in the XY data cubes. This problem was 
solved by identifying continuum sources observed in 
the raw image with sources catalogued in the NRAO 
VLA Sky Survey (Condon et al. 1998) and perform- 
ing astrometry on these sources. 

For each galaxy, data cubes were constructed us- 
ing both natural and robust weighting. We present 
in Figure |^ the data cubes in the form of channel 
maps derived by natural weighting, because natural 
weighting appeared to strike the best balance between 
the signal-to-noise ratio and resolution. This fig- 
ure presents the channels containing signal for each 
galaxy, bracketed by several noise channels. The 
channel velocity for each subframe is given in the 
upper left corner. For most of the targets, which 
had emission in a large number of channels, only ev- 
ery other channel is depicted. For UGC 6429 and 
NGC 4688, which were detected in a smaller number 
of channels, every channel containing signal is shown. 

Following this standard visibility reduction, the 
data image cubes were transferred for further reduc- 
tion in the GIPSY environment (van der Hulst et 
al. 1992, GIPSY 2000). Copies of the data cubes 
were then smoothed to provide easier separation of 
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regions with signal from regions with no signal. Data 
in these copies were cropped at the 3a level. The 
cropped regions were then mapped back to the un- 
smoothed cubes, which then consisted only of regions 
containing signal. Following correction for primary 
beam, these regions were integrated along the veloc- 
ity axis into zeroth, first, and second moments of the 
flux, thus creating three spatial maps, representing 
the H I flux, bulk motion, and velocity dispersion re- 
spectively. To eliminate spurious noise peaks during 
this integration, a "windowing" function was used, 
such that only peaks seen in two or more consecutive 
channels were included in the integration. The H I 
flux was then scaled to the column density N of the 
gas, assuming optically thin H I, via the relation 

/S 
— (1 + v/cfdv (1) 
ab 

where S is the total flux density in Jy/Beam, and a 
and b are the major and minor FWHM axes of the 
Gaussian beam in arcseconds. 

Following construction of the moment maps, simu- 
lated single-dish line profiles for each galaxy were ob- 
tained by integrating the data cubes over the spatial 
axes. Total integrated flux measurements were ob- 
tained by integrating the line profiles along the veloc- 
ity axis. They were then compared to single-dish data 
obtained from various sources. These single-dish data 
were corrected for pointing errors and beam smearing 
following the two-Gaussian distribution model of He- 
witt et al. (1983), with the Gaussian parameters de- 
rived from the synthesis data obtained here by fitting 
a double Gaussian to the moment zero maps. The 
comparison of synthesis to corrected single-dish flux 
is presented as part of Table |[ Following the discus- 
sion of Hewitt et ai, we expect between 15% and 20% 
uncertainty in the derived corrected flux. In general, 
there is agreement between the total flux measure- 
ments derived here and the single-dish data within 
the errors. 

The first order moment maps were used to obtain a 
kinematic model of the galaxies by iteratively apply- 
ing the standard GIPSY routine ROTCUR (Begeman 
1989). This routine divides the galaxy into a number 
of rings, and fits a function 

V{x,y) = Vsys + VrotCOs6smi (2) 

where 

Q ^ -jx- xq) sin(r) + iy- yo) cos(r) 
r 



to the velocity field in each ring, allowing the system- 
atic velocity Vsys, the rotation velocity Kot, inclina- 
tion i, position angle F, and kinematic center (xo^yo) 
to vary. 

Because of the face-on aspect of the target galax- 
ies, adequate kinematic solutions could be found, but 
were spread over a large area of parameter space. 
Specifically, although Vrot sin i remained well-constrained 
for each ring, changes in Vrot and i were found to 
be interchangeable with little effect on the quality of 
fit in the model function. The result was often that 
the algorithm produced unphysical values of Vrot with 
wildly varying inclinations if all parameters were al- 
lowed to vary independently. For this reason, the 
ROTCUR model was applied iteratively to the data, 
alternately holding Vrot{r), i{r), and {xQ,yQ) fixed. 
We present rotation curves and all subsequent analy- 
sis in terms of V^of sini. We caution, therefore, that 
variations in the VrotS\ni{r) curves we present may 
represent either actual changes in the rotation veloc- 
ities of the galaxies, changes in the galaxies' inclina- 
tions, or a combination of both. Inability to further 
constrain Vrotif) also precluded the construction of 
reliable model velocity fields which could reasonably 
be expected to reproduce the potentials of the galax- 
ies. Once the tilted ring fit was applied iteratively to 
the entire galaxy to obtain initial values of the pa- 
rameters, the model was fitted to both halves of the 
galaxy separately, so that Vrotsaii{r) and position 
angles could be determined for both sides indepen- 
dently. 

Position-velocity diagrams for one beam-width 
along the major axis of each data set were then ob- 
tained. To construct the position-velocity diagrams, 
slices of the data cubes were taken through the kine- 
matic centers obtained in the model fitting. Each slice 
was a strip of data whose spatial length was the total 
detected extent of the galaxy, spatial width was the 
width of a beam, and frequency depth was the depth 
of the cube. The spatial position angle of the slice 
was determined by averaging the model position an- 
gles over radii r < R25 . Data were then averaged over 
the width of the slice to produce the position-velocity 
diagrams. 

Distances to most of the sample galaxies were de- 
termined using a nonlinear Virgocentric infall model 
with an infall flow velocity of the Local Group of 
230 kms^^, Hubble parameter iJo — 74kms~^ Mpc~^, 
distance to Virgo of 16.6 Mpc, and the mean over- 
density of the Virgo cluster compared to the ambient 
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Table 3 
Global H I Properties 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Galaxy 


i^VLA 


FsD 




Vhi 


DH1/D25 


Vrot sin i 


Dist 


Mhi 


Mhi/Lb 




(Jy km 




(km 






(km s-i) 


(Mpc) 




(Mq/Lq) 


NGC 991 


21 


21^ 


91 


1531 


1.9 


33.0 


20 


2.0 


.27 


NGC 1042 


55 


56^ 


116 


1372 


1.7 


43.7 


20 


5.2 


.30 


UGC 3685 


43 


55^ 


97 


1796 


2.1 


36.8 


30 


9.1 


.54 


NGC 3596 


29 


33*^ 


127 


1193 


1.1 


51.7 


14 


1.3 


.20 


UGC 6429 


10 


11^= 


69 


3727 


1.5 


23.3 


75 


13 


.42 


NGC 4688 


28 


35^ 


60 


984 


1.8 


18.1 


17 


1.9 


.41 


NGC 5474 


78 


120^ 


60 


277 


1.9 


8.2 


7 


0.9 


.40 


NGC 5701 


58 


63'^ 


139 


1509 


2.6 


57.1 


26 


9.3 


.42 


UGC 12732 


79 


79^ 


139 


746 


3.1 


54.6 


15 


4.2 


3.2 



'^Fisher & TuUy 1981 
''HHMRvZ 

'^Staveley-Smith & Davies 1988 



density of pvirgo/po = 2. Our model uses these pa- 
rameters to solve for the velocity field of the local 
universe, then to find the possible distances of indi- 
vidual galaxies by finding the roots of the difference of 
the Local Group-ocentric velocity (derived from the 
heliocentric velocity by the method of Mould, Aaron- 
son, and Huchra (1980)) and the velocity found for a 
given distance. This is done by using the Numerical 
Recipes (Press et ai 1986) routine ZBRAK to deter- 
mine the bounds of the regions surrounding each root, 
which are then used to solve for the roots themselves. 

This model yields multiple distances for galaxies 
in the triple-valued-region close to the Virgo clus- 
ter. Distances of 12, 14, and 24 Mpc were ob- 
tained for NGC 3596, and of 7, 17, and 22 Mpc for 
NGC 4688. Here, we have taken the value of 14 Mpc 
for NGC 3596, based on the distance modulus of 31.41 
found by Bottinelli et ai (1985) using the TuUy- 
Fisher relation of H I line width and -B-band luminos- 
ity. We have taken a distance of 17 Mpc for NGC 4688 
based on a study of supernovae by Gaskell (1992) in 
which a distance of 17 Mpc yields a peak luminosity of 
SN 1966B commensurate with other type II-L super- 
novae. Following RDH, the distance to NGC 5474 of 



7 Mpc was taken directly from Sandage & Tammann 
(1974), where distances to galaxies in the MlOl sys- 
tem were found by measuring the angular sizes of H II 
regions. We have also adopted a value of 20 Mpc for 
NGC 991, equal to that of NGC 1042, because both 
galaxies are found by Garcia (1993) to be members 
LGG 71, a loose group of galaxies which includes the 
NGC 1052 group. 

The global H I properties derived from the anal- 
ysis are presented in Table ^. Table columns repre- 
sent: (1) The NGC or UGC designation of the target; 
(2) The integrated total flux observed at the VLA, 
i^VLA, in Jy km s~^; (3) The single dish flux, Fsd, in 
Jy km s^^, corrected for pointing errors and beam di- 
lution as described earlier in this section with associ- 
ated reference for the uncorrected flux measurement; 
(4) The full width of the H I line profile measured at 
20% of the peak flux, W20, in km s~^; (5) The system- 
atic velocity Vhi of the target, in km s^^, obtained 
from tilted ring modeling of the data, with typical 
errors of ±5 km s~^; (6) The ratio of the H I to opti- 
cal diameters, DH1/D25, where Dhi is measured at a 
surface density of IMq pc~^. (7) The observed value 
of Vrot sin i averaged over the flat part of the curve; 
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(8) The distance to the target, in Mpc, obtained from 
the Virgocentric infall model; (9) Total H I mass ob- 
served at the VLA, obtained from the flux density via 
the relation 

where D is the distance to the galaxy and 5*^ is the 
H I flux density in Janskys; and (10) the H I mass to 
blue luminosity ratio Mri /Lb, where the luminosity 
is derived from Bj,. 

2.3. Qualitative Features of the Target Ob- 
jects 

General comments regarding the gas morphology 
and dynamics of each target are discussed in this sub- 
section. The text describing each galaxy is accompa- 
nied by a figure presenting, in the top two frames, the 
H I column density contours overlaid on the optical 
image of the galaxy obtained by KHL, and a grayscale 
map of the H I column density itself. Noise in the col- 
umn density maps differs from pixel to pixel due to 
the windowing function, and was corrected for noise 
correlations among Banning smoothed channels fol- 
lowing the discussion of Verheijen (1997). The noise 
in a given pixel is 

= ^[[N -i/A) + N^N^Y <7^ (5) 
v6 

where N is the number of adjacent Hanning smoothed 
channels which contributed to the pixel, <7^ is the 
noise in a single channel, and 

is a factor representing the noise contributed by sub- 
tracting the first iVi and last N2 line-free channels 
from the cube. The average pixel value for all pixels 
with signal-to-noise of 2.75 < S/N < 3.25 was taken 
to define the 3a level for the map. Column density 
contours in the maps are selected for each galaxy in- 
dividually based primarily on clarity of presentation, 
and the H I grayscales are continuous and linear. 

In order to produce the optical overlays, astrom- 
etry was performed for each of the i?-band images 
obtained in KHL corresponding to the VLA target 
galaxies, based on star positions obtained from the 
HST Guide Star Catalog (Lasker et al. 1992). This 



also allowed morphological asymmetry measurements 
to be made with respect to the position of the galax- 
ies' optical nuclei. 

The center left frame exhibits the position-velocity 
diagram along the H I dynamical major axis, as de- 
termined by averaging over the position angles for all 
radii inside of i?25 . The position angle is given at the 
top of the plot, and the position axis is measured in 
arcminutes from the dynamical center. Dotted ver- 
tical lines in the position-velocity diagram indicate 
i?25. Graininess of the archive data is due to large 
channel widths and small number of channels in these 
sets, which make the corresponding diagrams data- 
sparse. 

The center left frame is a grayscale representa- 
tion of the H I velocity distribution, with grayscale 
steps equal to the contours. Again, the contours 
and grayscales are chosen individually for each galaxy 
based on clarity of presentation. Dark grayscales 
represent the receding side of the galaxy; lighter 
grayscales the approaching side. Values for selected 
contours are presented in km s^^. Overlaid on the 
grayscale moment frames are circles of radius i?25 
centered on the center of light of the optical galax- 
ies, except for NGG 5474, where the center of light is 
so displaced from both the dynamical center and the 
center of the outer optical isophotes that the dynam- 
ical center is used. We have elected to show circles 
here, although the RG3 indicates some elongation, as 
evident from Table |l|, because of ambiguity in extract- 
ing asymmetry versus inclination for face-on objects. 

Finally, the lower two frames present the integrated 
H I line profile and K-otsini(r) for the galaxy. This 
latter diagram presents the data for the line-of-sight 
rotation velocity Vrotsnii (km s~^), and position an- 
gle (degrees east of north) from the tilted ring fit- 
ting algorithm as a function of radius in arcseconds. 
The error values given are the numerical errors in the 
fitting algorithm; we estimate that systematic errors 
are on the order of twice as large. The data for the 
approaching side (circles; solid lines) and the reced- 
ing side (triangles; dashed lines) are presented super- 
posed. It is important to note that the ordinate of 
the rotation plot does not start at zero; because of 
beam smearing, we cannot trace the rotation close 
to the center. The differences between values for the 
approaching and receding sides would be washed out 
with an increased scale. 



9 



2.3.1. NGC 991 

NGC 991 (Figure |) is a flocculent SBc galaxy 
whose northwestern spiral arms extend twice as far 
from the nucleus as do the visible southeastern arms. 
Since the optical flux of the flocculent structures 
northwest of the nucleus is greater than that to the 
southeast, it is classified as asymmetric by KHL. The 
bar is elongated in the east-west direction. The de- 
tected neutral hydrogen extends to 1.5 times the opti- 
cal radius. A neutral hydrogen disk of 1.2 x 10^^ cm~^ 
coexists with the optical disk, except for a gap over 
the underluminous southern disk which has an aver- 
age H I column density of 8.5 x 10^*^ cm^^. This gap 
appears in a partial H I ring with two flux maxima at 
a radius of 1' (6.4 kpc). The outer H I disk is morpho- 
logically symmetric. According to Garcia (1993), the 
nearest neighbor to NGC 991 in LGG 71 is NGC 1022, 
34' (198 kpc at the adopted distance of 20 Mpc) to 
the northeast. 

The H I dynamics of NGC 991 are relatively undis- 
turbed within the optical radius as illustrated in the 
position-velocity diagram (see Figure ^ , but certain 
nonaxisymmetric dynamical features are evident in 
the outer disk. Closed contours on the east side of 
the major axis indicate a possible warp in the disk. 
Strong twisting of the minor axis velocity contours 
analogous to the "S-" shaped velocity contours of 
NGC 5474 (see § |2.3.7D is also present in NGC 991. 



The same phenomenon can also be seen in the chan- 
nel maps as an area of enhanced flux to the east of the 
galaxy at 1513.8 and 1508.6 km s~^. This is indicative 
of either a severely warped disk or strong streaming 
motions through an outer spiral arm. The dynamical 
center of the velocity distribution is consistent with 
the 7? -band center of light of the galaxy, but the dy- 
namics exhibit a 40° maximum shift in position angle 
from approaching to receding side for the same ra- 
dius, as a result of the detected "S-" shaped velocity 
contours. The synthesised line profile is asymmetric, 
with decreasing flux towards the approaching side. 

2.3.2. NGC 1042 

NGC 1042 was first mapped in H I by van Gorkom 
et al. (1986) whose primary purpose was to obtain 
H I synthesis data for the elliptical NGC 1052 and 
its companions, which include NGC 1042, some 15' 
(87 kpc at the adopted distance of 20 Mpc) south- 
west of NGC 1052. They point out that the neutral 
gas in NGC 1042 is asymmetric in the direction of 



NGC 1052, suggesting a tidal interaction. Because 
these data were collected with only 18 antennae yield- 
ing a synthesised beam of 60" x 60" at a velocity reso- 
lution (41 km s^^) more suited to the inclined galaxy 
NGC 1052, it was decided that newer data were re- 
quired to better constrain the dynamical parameters. 
These new data we present here. 

NGC 1042 (Figure ^ is a nearby asymmetric spiral 
galaxy whose eastern spiral arm contains many active 
H II regions, but whose western spiral arm is under- 
luminous. Although the total H I fiux agrees with 
previous observations, the H I mass derived from our 
data is larger by a factor of two than that reported 
by van Gorkom et al. and references therein, due to 
the difference in adopted distances. 

The H I morphology of this galaxy is unremarkable 
inside the optical radius, with the H I column density 
following the pattern of the dominant spiral arms out- 
side of the central H I hole. Asymmetry in the outer 
regions, however, is clear. The column density drop- 
off towards the northeast is far more gradual than 
the precipitous drop-ofl' to the southwest, giving the 
galaxy the appearance of compression on the south- 
western "forward" edge, as if traveling parallel to its 
disk through an intergalactic medium. Interestingly, 
this would indicate that NGC 1042 is traveling in 
the direction opposite to the position of NGC 1052. 
Van Gorkom et al. also comment that the H I in 
NGC 1042 is not significantly larger in extent than 
the optical galaxy. Our estimate of DH1/D25 — 1.7, 
however, is indicative of a large H I disk and is close 
to the average for our sample. 

The dynamics of NGC 1042 include substantial ev- 
idence along the minor axis of a warp in the disk, 
especially north of the optical disk. Several other 
nonaxisymmetric kinematic features are indicated as 
well. One of these is the closed velocity contours on 
the eastern side of the galaxy, seen in both the ve- 
locity field and the channel maps, implying a peak 
in Vrot just inside the optical radius. This feature is 
also consistent with a warp. The warped disk is also 
manifest in the H I position angle curve, in which T 
changes by as much as 43° as a function of radius. 
Streaming motions are apparent on both sides of the 
minor axis, correlating with the outer extensions of 
the two dominant optical spiral arms. Vrot siiii{r) in 
NGC 1042 is reminiscent of rotation curves presented 
by S3vA, in which opposite sides rise with different 
slopes. These dynamical perturbations are suggestive 
of tidal interaction with NGC 1052. 
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An interesting feature to note is that the position 
angle of the optical galaxy reported by KHL, found 
by fitting elliptical isophotes, appears perpendicular 
to the major axis of the velocity field of the H I. This 
effect is apparent in Figure |3| where the morpholog- 
ical major axis of the low-surface brightness optical 
galaxy appears to run north-south, while the dynam- 
ical major axis in the H I velocity field is clearly 
east-west. The implication is that the optical disk 
possesses inherent ellipticity, again suggesting tidal 
interaction with NGC 1052. Since KHL estimated 
inclinations based on optical axial ratios and the as- 
sumption of circular disks, the optical asymmetry pa- 
rameter derived there is likely to be an underestimate. 
In this work, we use the dynamical position angle. 

2.3.3. UGC 3685 

UGC 3685 (Figure |) is a flocculent SBb galaxy 
whose bar terminates at a well-defined circular ring. 
Two major spiral arms extend outwards from the ring, 
but their positions are uncorrelated with the bar. The 
northern spiral arm and associated flocculent struc- 
ture is more luminous than the southern structure, 
and contains many bright H II regions. An additional 
area of low surface brightness extends to the north- 
east. It is classified as symmetric by KHL. There are 
no galaxies within 1° and 1000 km s'^ of UGC 3685 
listed in the AGC. 

The neutral hydrogen distribution extends to a ra- 
dius of 2.1 i?25. Not all of the single-dish flux was re- 
covered during these observations, and although com- 
parison with an HI line profile obtained with the for- 
mer 91-m telescope at Green Bank available in the 
digital archive of H I spectra of Giovanelli & Haynes 
indicates that much of the difference can be explained 
by calibration errors, additional low-surface-density 
diffuse gas whose flux was not detected at the VLA 
cannot be ruled out. The "fishtail," a small western 
extension of the distribution, occurs at the 10^° cm~^ 
level, and along with the H I gap to the northwest and 
strong spiral structure, significantly affects the out- 
come of the morphological analysis. The overall dy- 
namics of the neutral hydrogen seem relatively undis- 
turbed, although streaming motions of 3-5 km s~^ 
are apparent just outside the optical radius to the 
north, and at 1.5i?25 to the south, following the spi- 
ral arms. The rotation curve and line profile of the 
galaxy, however, are asymmetric between approach- 
ing and receding sides. 

The most interesting dynamical feature of UGC 3685 



appears in the position-velocity diagram and the 
channel maps. A region of high (40 km s~^) veloc- 
ity dispersion is detected 2' northwest of the galactic 
nucleus at an azimuth of approximately 299°. This 
feature is consistent with the presence of a superbub- 
ble, an outflow of hot gas into the halo as a result of 
rapid star formation and supernovae type II in large 
H II regions. The symmetry in the velocity contours 
above and below the systematic rotational velocity of 
the remainder of the disk implies that the bubble orig- 
inated with an event near the midplane of the disk. 
A similar H I bubble has been observed in Ml 01 by 
Kamphuis et al. (1991). The bubble in UGC 3685 
is observed nearly normal to the disk, allowing the 
enclosed kinetic energy of expansion into the halo to 
be directly measured. The H I mass enclosed in the 
bubble is 1.22 x lO^M©, resulting in an approximate 
kinetic energy of 2 x lO^'^ ergs, of the same order as 
energies found in superbubbles in MlOl. 

2.3.4. NGC 3596 

NGC 3596 (Figure ||) is an optically asymmetric 
Sc galaxy which exhibits two "grand design" spiral 
arms surrounded by a low surface brightness system 
of flocculent structure. The contrast of the optical 
image in the accompanying figure was selected to em- 
phasize the fiocculent outer structure so that the rela- 
tionship between this structure and the H I morphol- 
ogy is clearer. The neutral gas in NGC 3596 is dis- 
tributed in a ring comparable in radius to the bright 
spiral arms, from which extends high-column density 
spiral structure in the south, while a single H I spi- 
ral arm extends to the north. The northern arm and 
the highest column density southern arm correspond 
to very low-surface brightness features in the optical; 
however, the optical features do not extend coherently 
along the entire length of the gas features. While the 
northern H I arm overlies its corresponding flocculent 
optical arm, the optical feature associated with the 
greatest H I column density in the southern arm is 
offset from the H I by 10" (0.7 kpc) to the east. The 
remaining H I spiral arms overlie corresponding low- 
surface brightness flocculent structures. The overall 
asymmetric distribution of neutral gas contributes to 
a relatively high value of the morphological asymme- 
try parameters for this target. 

Although there are no galaxies within 1° of NGC 3596 
with known redshifts near 1193 km listed in the 
AGC, NGC 3596 is dynamically perturbed. Although 
the H I disk does not extend far beyond the optical 
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radius, and the inner velocity contours of the galaxy 
appear very regular, warping is evident beyond R25, 
and the dynamical position angle of the approaching 
side varies southward by as much as 31 degrees. The 
severity of the warp is most evident in the channel 
maps, in which the disk appears so warped that sev- 
eral channels on both the approaching and receding 
sides contain flux in a circular or multiple-component 
pattern. Closed contours are evident on both the re- 
ceding and approaching sides, and the galaxy exhibits 
a falling Vrot sinj(r), indicating a warp. The line pro- 
file of the galaxy, however, is quite symmetric, hiding 
the inherent dynamic nonaxisymmetry, in agreement 
with HHMRvZ. 

2.3.5. UGC 6429 

As measured by KHL, UGC 6429 (Figure |) is opti- 
cally an asymmetric Sc disk exhibiting two inner and 
two outer spiral arms, with the southern arm termi- 
nating in a large, bright H II region. The H I dis- 
tribution exhibits an asymmetric ring following the 
inner spiral arms, and a large high-column density 
feature to the east, just outside R25. This feature is 
an extension of the northern optical outer spiral arm, 
which gives the entire distribution a very asymmetric 
qualitative appearance. 

Dynamically, the most apparent departures from 
axisymmetry in UGC 6429 are the closed contours 
on the approaching side of the velocity field, and the 
change in position angle with radius to the north. 
There is no dynamical feature corresponding to the 
eastern morphological feature. The asymmetry in 
the integrated line profile of this galaxy is due to 
the nonuniform weights given in integrating the line 
profile due to the eastern density enhancement. The 
nearest neighbor listed in the AGC to UGC 6429 is 
NGC 3668, 17' (371 kpc in projected distance) to the 
southeast, at 3508 km s^^. 

2.3.6. NGC 4688 

NGC 4688 (Figure ^ is an outlying member of the 
Virgo cluster, and is a member of LGG 292 (Gar- 
cia 1993). Its nearest neighbor listed in that group 
is UGC 7983, an irregular galaxy 43' (213 kpc at 
the adopted distance of 17 Mpc) to the southeast. 
NGC 4688 is a very asymmetric Scd galaxy in the op- 
tical, where flocculent spiral arms extend away from 
the main body of the galaxy to the northeast. Besides 
the tidally deformed galaxy NGC 5474, NGC 4688 is 



the most optically asymmetric galaxy in the sample. 
The H I morphology, too, exhibits strong asymme- 
tries, particularly the very high-column density re- 
gions to the east-southeast just at R25 . These regions 
correspond to the southernmost of the flocculent opti- 
cal arms as seen in i?-band. Spiral structure, in fact, 
can be seen throughout the H I disk, corresponding to 
flocculent structures in the optical galaxy. The low- 
column density H I is distributed preferentially to the 
northeast, lending an even more qualitative lopsided 
appearance to the galaxy. This feature contributes 
to a larger value of RH1/R25 to the north than to 
the south. Quantitatively, as calculated in the next 
section, the morphological asymmetry of the H I dis- 
tribution makes NGC 4688 the most morphologically 
asymmetric galaxy in gas as well as in optical light. 

Dynamically, too, NGC 4688 is disturbed. The 
receding side of the velocity field exhibits closed con- 
tours which reopen again at 1.5i?25, a feature simi- 
lar to that found in NGC 991 and NGC 5474. This 
feature leads to a large asymmetry in the rotation 
curve, which manifests itself as a large difference be- 
tween observed velocity on the approaching and re- 
ceding sides at the same radius. The global profile is 
unusual too, in that it exhibits very gradual declines 
from the peaks to the noise level on both sides of the 
profile. 

A possible explanation for this very asymmetric be- 
havior lies 6.'7 (33 kpc) to the northeast of the dynam- 
ical center of NGC 4688: a smaU 3.2 x lO^M© com- 
panion detected in the H I map. The peak of the H I 
distribution of the companion is located at B1950.0 
coordinates a = 12''45'"26!5, S = +04°42'10". The 
systemic velocity of the peak of the companion's emis- 
sion is 991 km s~^. 

The column density map of the companion is pre- 
sented in Figure ||. The top left frame of this fig- 
ure presents the column density of NGC 4688 and 
the companion for scaling; the top and bottom right 
frames detail the companion itself. The H I column 
density contours in the top right panel are overlaid on 
the Palomar Observatory Sky Survey]] plate for this 
region digitized by the Space Telescope Science Insti- 
tute Digitized Sky Survey^ The lower right panel 

^Tho National Geographic Society - Palomar Observatory Sky 
Atlas (POSS-I) was made by the California Institute of Tech- 
nology with grants from the National Geographic Society. 

*The Digitized Sky Surveys were produced at the Space Tele- 
scope Science Institute under U.S. Government grant NAG W— 
2166. The images of these surveys are based on photographic 
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depicts the POSS plate alone. The companion is just 
resolved by our observations in the CS array, and the 
maximum column density observed is 6 x 10^*^ cm~^. 
An optical counterpart to the H I emission is clearly 
seen. The observed location of the H I emission, how- 
ever, is offset 20" to the southwest of the maximum of 
the optical emission, but does correspond to a smaller 
local maximum. We identify these two optical max- 
ima as the galaxy pair CGCG 043-029. The bottom 
left frame of Figure ^ illustrates the synthesised line 
profile of this object. The maximum flux occurs in 
the 991 km s^^ channel. The gradual slope on the 
receding side of the profile, however, suggests that 
there is probably significant additional flux beyond 
our velocity limit of 1040 km s^^. 

An H I detection from the northeastern component 
of CGCG 043-029 appears in a table of H I obser- 
vations of faint galaxies conducted with the Arecibo 
circular feed compiled by Lu et al. (1993), ex- 
cept that the heliocentric velocity reported there is 
4889 km s~^. It is a possibility, therefore, that the 
pair is an optical pair only, with the 4889 km s~^ emis- 
sion associated with the northeastern component, and 
the 991 km s^^ emission associated with the south- 
western component. The data, however, are too in- 
complete for us to make any supportable claims re- 
garding this object beyond its mere detection. 

During the initial analysis of the 21 cm line data, 
we imaged the data cube using a tapered beam of 
40 X 45 arcseconds size to examine the data for ev- 
idence of an H I bridge between NGC 4688 and 
CGCG 043-029. No such bridge was detected. 

2.3.7. NGC 5474 

NGC 5474 (Figure ||) is a peculiar spiral galaxy 
approximately 44' (90 kpc at the adopted distance of 
7 Mpc) to the south of MlOl, and is tidally deformed 
into a very asymmetric and disturbed morphology as 
observed by KHL in the optical. 

NGC 5474 was first observed in H I synthesis by 
van der Hulst & Huchtmeier (1979) at the WSRT, 
with significant H I fiux appearing in four channels 
with a velocity resolution of 27.2 km s~^. Subse- 
quently, this galaxy was observed by RDH with the 
D array at the VLA. Both groups indicate that the 



kinematic major axis of this galaxy does not appear 
to form a straight line, but rather is contorted into 
an "S-" shape, and also note that the kinematic cen- 
ter and morphological centroid of the H I distribution 
are approximately co-located with the center of the 
outer optical isophotes, and not with the offset opti- 
cal nucleus. Hence, the true center of the galaxy is 
not the bright "pseudonuclear" region located to the 
north of the overall optical light distribution. 

As mentioned in §2.1.2, we are able to trace H I 



data obtained using the Oschin Schmidt Telescope on Palomar 
Mountain and the UK Schmidt Telescope. The plates were 
processed into the present compressed digital form with the 
permission of these institutions. 



to a radius of approximately 400", in contrast to the 
250" radius recovered by RDH and by van der Hulst 
& Huchtmeier. In addition to the "clumps" of neutral 
hydrogen found in the high HI column density regions 
by RDH, we find a smooth HI distribution beyond 
twice the optical radius. Even so, only half of the re- 
ported single-dish flux for this galaxy was recovered, 
because NGC 5474 and MlOl share a common H I 
distribution (Huchtmeier & Witzel 1979), and even 
the short spacings of the D array are unable to re- 
cover the diffuse gas which envelopes this galaxy and 
its dominant companion. We do, however, detect 50% 
more flux than that found by RDH. 

Dynamically, the galaxy exhibits a number of bizarre 
features. RDH detect 50° of position angle change 
from the dynamical center to a radius of 250". Our 
reanalysis agrees with their analysis and also reveals 
the continuation of this trend further out, where the 
position angle varies by a total of 95° from the dy- 
namical center to a 400" radius. Position angles for 
the approaching and receding sides also diverge in the 
outer regions, with the difference in position angle 
from one side to the other of « 50°, although the 
greatest divergence occurs in a region of small signal- 
to-noise (slightly below 3a) and may therefore be an 
artifact of the model. These large position angle vari- 
ations as a function of radius are largely due to the 
very unusual dynamic feature found 40" to the south 
of the optical center of the galaxy disk, in which the 
velocity contours fully reverse their directions, a con- 
tinuation of the "S-" shape reported by the previous 
observers. This feature also displaces the dynamical 
center of the outer regions about 1' to the southwest 
of the optical pseudonucleus center of light. The fea- 
ture contributes to the measurable asymmetry in the 
integrated single-dish line profile. Such a feature may 
be the result a very warped disk viewed in a direction 
which bisects the angle formed by the fiat and warped 
portions of the disk. The pseudonucleus itself has no 
discernible effect on the dynamics of the disk. 
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RDH use a model in which the galaxy inclination 
is fixed at 21° to determine their rotation curve; if 
we fix the inclination at this value, we recover the 
same curve. Both analyses agree on the general shape 
of the curve, which rises out to about 100", steeper 
on the receding side, and flattens out with a shallow 
minimum at 200". 

2.3.8. NGC 5701 

The barred Sa galaxy NGC 5701 (Figure |l|) ex- 
hibits prominent inner and outer ring structures in 
the optical image along with a large bulge and bar at 
position angle —5°. The optical galaxy is embedded 
in a very large H I envelope of DH1/D25 — 2.6. More 
H I flux is observed on the west side of the galaxy 
than on the east by a factor of two, and the western 
drop-off is more precipitous than the east. This fea- 
ture is similar to that found in NGC 1042 (§ p.3.2D , 
where it is interpreted as a possible remnant of tidal 
interaction. In this case, however, no companion to 
NGC 5701 is seen in our optical or H I data. It is, 
however, identified by Garcia (1993) as a member of 
the group LGG 386, the nearest member of which to 
NGC 5701 is UGC 9385, 57' (431 kpc at the adopted 
distance of 26 Mpc) to the west. 

In this connection, it is important to note that 
NGC 5701 was also mapped at Arecibo by DuPrie 
& Schneider (1996). We agree with their Figure 3, 
that the H I is extended to ^ 10'; however, we do not 
recover the H I "cloud" they report 15' northwest of 
the center of NGC 5701 and of flux 5.5 Jy km s~^. 
We do detect two flux peaks 4(t above the noise 13' to 
the northwest in two channels (1495 and 1505 km s~^) 
in the archive data. However, because the emission 
peaks occur in two different spatial locations sepa- 
rated by about a beam-width, these peaks were ex- 
cluded from the moment maps by the windowing func- 
tion. That these peaks occur at the same heliocentric 
velocity as the "cloud" reported by DuPrie & Schnei- 
der, and would represent a total flux of 4.6 Jy km s~^ 
suggests that the "cloud" is marginally detected in 
the synthesis data. 

Nonlinearity of the velocity contours along the mi- 
nor axis of NGC 5701 can be observed in the veloc- 
ity field. This results in a dynamical position angle 
change of 43° as a function of radius, similar to the 
value obtained from analysis of NGC 1042. NGC 5701 
also exhibits significant differences between the val- 
ues of Vrot sini(r) on either branch of the major axis, 
again similar to the rotation curve of NGC 1042. 



NGC 5701 differs, however, in that it exhibits a deep 
trough in Vrotsmi{r) at a radius of 210". This is 
suggestive of either a warp or strong streaming mo- 
tions in this galaxy. The synthesised line profile of 
NGC 5701 is consistent with the single-dish profile 
presented by HHMRvZ. In agreement with them, we 
find that the profile is symmetric. 

2.3.9. UGC 12732 

UGC 12732 (Figure |l|), first observed by Schul- 
man et al. (1997), is a small Sm galaxy with a 
roughly symmetric optical appearance. Compared to 
the optical galaxy, the neutral gas structure detected 
is extremely large: DH1/D25 = 3.1 as given in Ta- 
ble ^, and Mhi /L b is an order of magnitude greater 
than that of the other targets. Indeed, UGC 12732 
is more gas-rich (or star-poor) than most Sm galax- 
ies; Mh/Lb of 3.2 is well in excess of the median 
of 0.78 for Sm and Im galaxies given by Roberts 
& Haynes (1994). The nearest neighbor of known 
redshift within 1000 km s'^ to that of UGC 12732 
is NGC 7741, 45' (196 kpc) to the southeast, at 
744 km s-i. 

The contours in Figure |ll| were chosen to empha- 
size the spiral structure slightly to the north of the 
optical galaxy. This spiral structure does not appear 
to be associated with any optical spiral structure. The 
region of maximum H I column density is located ap- 
proximately 2' to the SSE of the optical center, which 
corresponds to the termination of the eastern spiral 
arm as seen in the i?-band. The outer H I density 
contours are centered on the optical center of light, 
as is a 1' radius ring surrounding the inner regions of 
the optical galaxy. The spiral structures extend from 
this ring and are especially apparent to the north- 
west. These structures contribute to a very high value 
of morphological asymmetry for the neutral hydro- 
gen distribution. In dynamics, a 30° position angle 
change on both the approaching and receding sides of 
the galaxy, rising Vrot sin j(r), and an asymmetric line 
profile contribute to overall dynamical nonaxisymme- 
try in UGC 12732. 

As Schulman et al. commented, the large DH1/D25 
ratio may be an indicator of inefficient star formation 
in this galaxy. The Vrot sin i{r) curve and position 
angle fits we derive with the tilted ring model dif- 
fer from those derived by Schulman et al., in that 
the qualitative shapes of the curves significantly dif- 
fer between the two analyses. While in the Schulman 
et al. analysis, the inclination was allowed to vary to 
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obtain a flatter rotation curve, our analysis does not 
attempt to disentangle i from K-ot; resulting in a ris- 
ing Vrot sin «(r) curve and a greater change in position 
angle as a function of radius. 

3. Comparing Morphological and Kinematic 
Asymmetry 

3.1. Obtaining the Various Measures of Asym- 
metry in a GalELxy 

We examined our data sets for any correlations 
between asymmetric morphology and nonaxisymmet- 
ric dynamics. In doing so, we wished to search for 
methods of quantifying deviations from kinematic 
axisymmetry, including both non-planar and non- 
circular motions, which coidd then be compared to 
the morphological measures. Dynamic models of the 
galaxies allowed us to define several kinematic quan- 
tities which, we hypothesised, might be correlated 
with asymmetric morphology. The features we exam- 
ined were: morphological asymmetry in i?-band and 
H I, change of kinematic position angle with radius, 
change of kinematic position angle with azimuth, and 
asymmetry in rotation curves. Single-dish H I line 
profile asymmetry was also computed for these galax- 
ies using the method described by HHMRvZ. 

The zeroth order H I moment maps were analyzed 
for morphological asymmetry using the "Method of 
Sectors" described in KHL. Each galaxy was divided 
into 8 equal-area sectors with the center of the pat- 
tern located at the optical center of light, and the po- 
sition angle of the sector pattern oriented along the 
dynamical position angle of the galaxy. The total H I 
flux was measured in each of the sectors. The max- 
imum flux differences were then normalized by the 
total integrated H I flux through the pattern. The 
result is the asymmetry parameter Af. 

/\ Jmax 

where A/"^"^ is the largest flux difference between 
any two sectors, and /„ is the flux in the nth sector. 

ii-band values of A j^r are presented here adapted 
from the data used by KHL. In KHL, sector analysis 
of optical data was presented by taking ratios of the 
maximum flux in a sector to the minimum flux in a 
sector, and converting this to a magnitude difference 
in the standard way. Thus, KHL present the data 
in terms of magnitude differences AM between sec- 
tors rather than normalized flux difl[erences. There, a 



galaxy was said to be asymmetric if a value of AAf 
differed from zero by more than the ha RMS noise 
level, regardless of the magnitude of AM itself. For 
this reason, some galaxies with large values of AM 
were classified as symmetric because the errors were 
also large. Similarly, small values of AM with small 
errors could indicate asymmetric galaxies. Of the nine 
galaxies considered here, UGC 3685, NGC 5701, and 
UGC 12732 were classified as symmetric; all others 
were classified as asymmetric. Here, we use the same 
fiux data as were used by KHL, but present them 
as normalized flux differences Af for both the optical 
and H I data because they arc more useful in this form 
for understanding the symmetry properties of galaxy 
morphology in an absolute and directly comparable 
sense. 

Deviations from axisymmetric dynamics of the tar- 
gets was measured using data obtained from the 
model fitting. Because a symmetric, flat rotating 
disk should exhibit a constant position angle, we took 
changes in position angle in the model as quantifiers of 
nonaxisymmetry. The first of these is the maximum 
change in position angle on one side of the galaxy 
as a function of radius, Ar(r). Note that in cases 
where the position angle changes steadily with radius, 
Ar(r) can depend on the maximum radius for which 
we can construct a model. As applied to these face- 
on galaxies, Ar(r) is most sensitive to non-planar, as 
opposed to non-circular motions. The second of these 
is the diff'erence between position angles on approach- 
ing and receding sides of the model at the same radius 
averaged over all radii, Ar(^). Ar(^) is sensitive to 
non-circular motions and asymmetric warps. 

Extending the analysis of S3vA, asymmetry in 
Vrot sin i{r) curves was quantified in the following way. 
The curve was folded over the kinematic center of the 
model. The total area between the two curves was 
then measured and normalized by the average area 
under the curve, to produce the measure of global 
rotation asymmetry 82'- 

_ 2 / \\Vrot{r)smi{r)\ - |Kot(-r) sini(-r)|| dr 
J V2 [\Vrot{r)smi{r)\ + \Vrot{-r)smi{-r)\] dr 

(8) 

Clearly, for more inclined galaxies where Vrot and i 
can be fit separately, 5*2 c;an be modified to measure 
asymmetry in Vrot{r) alone. In our sample of face-on 
galaxies, 52 is very sensitive to non-planar motions. 
In more inclined samples, ^2 would be more sensitive 
to non-circular motions. 
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Because several groups (for example, Tifft & Cocke 
1988, Richter & Sancisi 1994, and HHMRvZ) have 
studied the symmetry properties of 21-cm line profiles 
as an indicator of dynamical asymmetry in a galaxy, 
here too we apply the method of line profile asym- 
metry described by HHMRvZ. The quantity defined 
there is: 



A 



SdV 



l/h 



(9) 



where Vi and Vh are the low and high velocity limits, 
respectively, of signal detection above 3cr, and Vmed 
is the median between them. The characteristic lop- 
sidedness of the line profile is given simply by 



A 



l/h, 



if Ai/h > 1 



1/A;//,, otherwise. 



(10) 



The values of A^ of 1.00 and 1.03 obtained in this 
work for NGC 3596 and NGC 5701, respectively, are 
in agreement with the values of 1/0.99 = 1.01 ob- 
tained by HHMRvZ within errors and differences in 
velocity resolution. HHMRvZ found that approx- 
imately 50% of the line profiles studied exhibited 
asymmetry parameters An > 1.05, which they con- 
sidered to be indicative of overall asymmetry. 

Results for the asymmetry measures for each tar- 
get are presented in Table ^. A plot of Af^fi versus 
Af^Hi and the plots of all morphological/dynamical 
pairs are presented in Figure |l^. This figure excludes 
points for NGC 5474, because of its large lower bound 
on Af^n. Errors on An are the size of the points, 
while systematic errors on S2 we estimate at about 
20%. Little correlation is apparent in our data set 
between H I morphology and dynamics for disk galax- 
ies in the field. In our sample, there is a relationship 
between the H I and i?-band morphologies in which 
large values of A/,_r imply large values of Af^Hi, al- 
though not the converse. It should not be surprising 
to find that optical and H I morphology are related, 
because H I spiral structure tends to follow the optical 
spiral structure. 

All pairs of morphological versus dynamical pa- 
rameters were also examined. The only statistically 
significant relationship between Af and any of the 
calculated dynamical parameters is a weak relation 
with Ar{9). Large values of AT{d) imply large val- 
ues of Af^Hi within our sample, but again, not the 
converse. It is also interesting that some of the most 
morphologically symmetric galaxies exhibit asymme- 
try in their line profiles as measured by An > 1.05. 



This qualitative observation and the lack of correla- 
tion between these parameters seems to support the 
casual observations of HHMRvZ and others who have 
noted a similar non-correlation. For instance, HHM- 
RvZ comment that very morphologically disturbed 
galaxies such as NGC 1637 often have undisturbed 
H I line profiles. 

The morphology of the asymmetries themselves is 
also worthy of study. KHL classified asymmetric mor- 
phology into three categories according to the relative 
locations of the maximum and minimum flux sectors 
in the 8-segment i?-band analysis. The categories 
were:: "bisymmetric," in which the extreme sectors 
were adjacent, "boxy," in which the extreme sectors 
were 90° apart, and "lopsided," in which the extreme 
sectors were 135° or 180° apart. Of the galaxies in this 
study classified as asymmetric by KHL, NGC 991, 
UGC 6429, and NGC 5474 were classified as lop- 
sided, while NGC 1042, NGC 3596, and NGC 4688 
were classified as boxy. We find that in the H I, all 
of these galaxies exhibit lopsided asymmetry, except 
for NGC 1042. NGC 1042 exhibits boxy asymmetry, 
caused by the peculiar perpendicularity of its mor- 
phological and dynamical major axes. The remaining 
galaxies in the current sample were classified as sym- 
metric by KHL, and therefore were not categorized in 
this way in i?-band. Of these galaxies, NGC 5701 and 
UGC 3685 are boxy in the H I, while UGC 12732 is 
lopsided. None of the observed galaxies were "bisym- 
metric" in the H L 

It is also of interest to search for asymmetry trends 
as a function of global gas parameters. Figure |l^ il- 
lustrates the relationships between Af^Hi and and the 
relevant global H I parameters of the targets as given 
in Table |^. Again, no clear relationships are appar- 
ent. The relationship between dynamical asymmetry 
and global H I properties is even more tenuous. Fig- 
ure ^ illustrates this clearly. There are no detectable 
dependencies of any global property of the gas on the 
rotation asymmetry 82- 

3.2. Comparison of Current Techniques 

As discussed in §|l|, several studies of optical mor- 
phological asymmetry have been conducted, in gen- 
eral using a number of different quantitative mea- 
sures. One of these, Conselice et al. (1999), includes 
NGC 3596 and NGC 5701 in its galaxy sample. They 
employ a method in which a galaxy image is rotated 
through some angle (j) and then subtracted from the 
original, unrotated image. The normalized residuals 
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Table 4 
Target Symmetry Properties 



(1) (2) (3) (4) (5) (6) (7) 

Galaxy Af^R Af^Hi AT{r) AT{e) S2 

(degrees) (degrees) 



NGC 991 


0.0254 ± 


.005 


0.0393 ± 0.002 


11.9 ± 2.5 


5.0 ± 


0.7 


9.9 


1.12 


NGC 1042 


0.0259 ± 


.005 


0.0260 ± 0.001 


43.4 ± 3.9 


5.1 ± 


0.4 


19 


1.04 


UGC 3685 


0.0536 ± 


.02 


0.0646 ± 0.001 


18.6 ± 3.9 


6.2 ± 


0.3 


15 


1.09 


NGC 3596 


0.0212 ± 


.001 


0.0491 ± 0.002 


27.2 ± 4.0 


6.6 ± 


0.4 


3.2 


1.00 


UGC 6429 


0.0189 ± 


.003 


0.0607 ± 0.003 


26.1 ± 5.0 


6.4 ± 


1.5 


9.4 


1.08 


NGC 4688 


0.0725 ± 


.005 


0.0690 ± 0.003 


25.9 ± 3.2 


9.0 ± 


0.8 


15 


1.07 


NGC 5474 


> 0.1 




0.0576 ± 0.01 


95.4 ± 6.4 


14.9 ± 


0.7 


10 


1.18 


NGC 5701 


0.0322 ± 


.01 


0.0419 ± 0.008 


43.0 ± 3.0 


4.1 ± 


0.3 


3.6 


1.03 


UGC 12732 


0.0297 ± 


.02 


0.()()l:-') ± 0.009 


1.").:-') ±0.7 


•1.5 ± 


1.4 


(i.O 


1.28 



are taken as the asymmetry measurement. Our mea- 
sure of R- band morphology is in disagreement with 
Consclicc et ai, insofar as NGC 3596 is one of our 
most symmetric galaxies as measured by Af^R, and 
NGC 5701 is measured to have moderate asymmetry 
by the same quantifier. In contrast. Consclicc et al. 
find that NGC 3596 is more asymmetric than about 
75% of their sample, and NGC 5701 is one of the most 
symmetric. These differences arc likely due to several 
factors; primarily, our methods are fundamentally dif- 
ferent, and measure different types of asymmetry. In 
addition, centering issues arc very likely to be an im- 
portant factor in the disagreement here, as Conselice 
et al. use a different centering method from that used 
here (see below). 

The dynamical symmetry properties which we have 
derived represent a wide range of dynamical variables, 
such that the combination of these parameters allows 
for the study of the global properties of a galaxy which 
are the most easily accessible to the casual observer. 
Other techniques for the determination of dynamic 
asymmetry have also been proposed, such as that 
by SFdZ. These authors present methods by which 
lopsided spiral modes and global elongation in the 
galactic potential may be inferred from the presence 
of nonaxisymmetric motions in the observed velocity 
field. In light of theoretical models such as that of Jog 
(1997), it is clear that perturbations in the gravita- 
tional potential can manifest as perturbations in the 



observed velocity fields. 

Unlike the work presented here, determining the 
spiral structure of the potentials of disk galaxies us- 
ing the method of SFdZ leads immediately to infor- 
mation about the total mass distribution in galactic 
halos. Nevertheless, the method is complex and is 
most applicable over a narrow range of inclinations. 
The parameters calculated in the present work are 
calculable over a wide range of galaxy inclinations 
and are quickly attainable from basic data, and thus 
are better suited to quickly identifying asymmetric 
galaxies from large samples. Once identified, more 
detailed methods would be invaluable for deriving the 
perturbed galactic potential of individual galaxies. 

A method based on SFdZ was used by S3vA to 
quantify the asymmetry in the two galaxies presented 
in that paper, NGC 4395 and DDO 9. There, they 
find that both galaxies, selected for strong kinematic 
lopsidcdncss, exhibit asymmetric rotation curves. To 
test this observation, we attempted to apply the SFdZ 
method to our sample of nine galaxies to find the 
zeroth through third harmonic components of each 
galaxy's velocity field as a function of radius. How- 
ever, because of the low inclinations of these galaxies, 
satisfactory decomposition into harmonic kinematic 
components for the comparison was impossible, as 
measured in part by consistently large values of the 
third cosine coefficient, 63, which should be zero when 
the inclination of the tilted ring fit is correct. 
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Closely connected with both the present method 
and those of Conselice et al., SFdZ and others is the 
question of centering. Although not yet utilized as a 
parameter of asymmetry itself, the question of find- 
ing the kinematic center and/or the center of light of 
galaxy is critical to the question of galactic symmetry; 
imprecise determination of the center of the galaxy 
can easily introduce unacceptable errors in the de- 
termination of symmetry properties. In the present 
method and that of Conselice et al, improper cen- 
tering creates a false dipole in the overlay pattern 
(or the rotation) for determining morphology, as it 
does for true Fourier morphology techniques such as 
that of Zaritsky & Rix (1997). Improper centering 
also contaminates measures of 52 and AF. In this 
work, centers for kinematic parameters were chosen 
iteratively for best chi-square fit in the tilted ring 
model, while the centers for the morphological pa- 
rameters were chosen as the centers of light of the 
optical galaxy. Conselice et al. iteratively find the 
centering which minimizes their morphological asym- 
metry parameter. SFdZ recommend selecting centers 
which minimize the total residual velocity. 

Beauvais & Bothun (1999) also consider the is- 
sue of centering when applying a tilted-ring model to 
Fabry-Perot data. They use three methods for cen- 
tering of the galaxy: photometric ellipse fitting, find- 
ing the center of symmetry of the tilted ring model, 
and fitting the center as a free parameter in the tilted- 
ring model. They report that each of these centering 
methods yielded equally good model velocity fields, 
but that it was common for the center to vary consid- 
erably in the outer radii of the galaxy when allowed to 
vary freely, and interpret this as possibly indicating 
actual misalignment of the inner and outer rotation 
axes. 

Because centering is so critical, nearly all studies 
of asymmetry must choose one or two methods for 
finding the center of a galaxy. These derived centers 
themselves can also be used as indicators of lopsid- 
edncss. Offsets of the kinematic center from the op- 
tical center of light and dependence on the derived 
kinematic center on radius can both be indicators of 
lopsided dynamics. In this work, the sample galaxies 
were chosen on the basis of their face-on inclinations 
for the purpose of examining the optical morphology, 
making the kinematic centers not sufficiently certain 
for comparison with morphological centers. However, 
future work with more inclined galaxies may be able 
to use centering itself as a symmetry indicator. 



4. Summary 

This work represents a concurrent analysis of both 
optical and 21-cm symmetry properties of a sample of 
disk galaxies. We have presented methods by which 
the symmetry properties of a galaxy can be quickly 
computed over a wide range of dynamic variables. 
These parameters are sensitive to both non-circular 
and non-planar motions, depending in large part on 
the inclinations of the galaxies. In this case, the face- 
on inclinations of the targets make the parameters 
most sensitive to non-planar motions. 

In our sample presented here, tidal interaction 
with companion objects dominates the asymmetry in 
NGC 5474 and cannot be ruled out for NGC 1042 
and NGC 4688. Nevertheless, dynamical nonaxisym- 
metries clearly prosper in more isolated and even mor- 
phologically symmetric galaxies such as UGC 12732 
and UGC 3685 as well. Pervasive deviations from 
axisymmetry even outside of tidal interactions seem 
to indicate that asymmetric galaxies are stable and 
long-lived phenomena. Our data seem, however, to 
indicate that there is little or no connection between 
lopsided morphology and non-planar motions in such 
galaxies. This lack of clear connection between mor- 
phology and warping seems to imply that separate 
physical mechanisms may be at work in each domain 
to produce deviations from flat, axisymmetric disks. 

The methods presented in this paper can be quickly 
applied to large samples of galaxies to identify devi- 
ations from axisymmetry in dynamics and morphol- 
ogy. The small sample size of this initial study and 
the kinematic uncertainties produced by the sample 
galaxies' face-on inclinations can only rule out the 
strongest of correlations among the various parame- 
ters. Larger samples and samples at varying inclina- 
tions are required to truly study the relationships be- 
tween these parameters at more than a cursory level. 
A future study with galaxies of a wider range of incli- 
nations will also allow for additional parameters such 
as kinematic offsets which can then be adequately 
studied. Future numerical work can also benefit from 
an easily applied method to track and identify de- 
partures from symmetry. Such numerical work could, 
for instance, study the origins and longevity of asym- 
metries, and predict the course of their evolution in 
spiral galaxies. 
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Fig. 1. — Channel maps for each of the six new data 
sets in the sample: (a) NGC 991; (b) NGC 1042; 
(c) UGC 3685; (d) NGC 3596; (e) UGC 6429; 
(f) NGC 3688. Channel velocity is given in each sub- 
frame, and contours are at the (-2, 2, 3, 6, 9, . . .)a 
levels, where a is as given in Table Grayscales are 
at the (6, 12, 18, . . .)a levels. 



Fig. 2.— VLA H I synthesis results for NGC 991. Top 
left panel: H I column density contours detected at 
the VLA, with light gray contour at 3.4 x 10^^ cm~^ 
and solid contours at {10.2,23.8,37.4,51.0,...} x 
10^^ cm~^, overlaid on the iZ-band image obtained by 
KHL. The first solid contour represents the Scr noise 
level. Top right panel: The same H I column den- 
sity contours as in the top left panel, plotted over a 
linear grayscale column density map. The circle is of 
diameter D25 and is centered on the iZ-band center of 
light of the galaxy. Center left panel: The position- 
velocity diagram obtained from the H I data cube. 
The slice was made along the dynamical major axis 
r = 295° determined by averaging the dynamical po- 
sition angle over radii r < i?25- The x-axis indicates 
arcminutes northwest of the dynamical center, and 
the dotted vertical lines identify i?25. Center right 
panel: The H I observed velocity field with contours 
and grayscales from 1490 km s~^ (east, light shading) 
to 1570 km (west, dark shading) in increments of 
10 km s-^ The 1500, 1540, and 1580 km s^^ con- 
tours are identified in the panel. The circle is of di- 
ameter D25 and is centered on the i?-band center of 
light of the galaxy. Bottom left panel: The synthe- 
sized H I line profile obtained by integrating the data 
cube along both spatial axes. Bottom right panel: 
The rotation Vrot sin i and dynamical major axis po- 
sition angle F as a function of radius in arcseconds, 
obtained from the tilted-ring model. The data for the 
approaching side (circles; solid lines) and the reced- 
ing side (triangles; dashed lines) are presented super- 
posed. 
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Fig. 3.— VLA H I synthesis results for NGC 1042. 

Top left panel: H I column density contours detected 
at the VLA, with light gray contour at 2.6 x 10^^ cm~^ 
and solid contours at {7.8,23.4,39.0,54.6,...} x 
10"'^'' cm^^, overlaid on the R band image obtained by 
KHL. The first solid contour represents the 3a noise 
level. Top right panel: The same H I column den- 
sity contours as in the top left panel, plotted over a 
linear grayscale column density map. The circle is of 
diameter D25 and is centered on the i?-band center of 
light of the galaxy. Center left panel: The position- 
velocity diagram obtained from the H I data cube. 
The slice was made along the dynamical major axis 
r = 292° determined by averaging the dynamical po- 
sition angle over radii r < i?25- The cc-axis indicates 
arcminutes northwest of the dynamical center, and 
the dotted vertical lines identify i?25- Center right 
panel: The H I observed velocity field with contours 
and grayscales from 1320 km s~^ (east, light shading) 
to 1420 km s""'^ (west, dark shading) in increments of 
10 km s-^ The 1330, 1370, and 1410 km s'^ con- 
tours are identified in the panel. The circle is of di- 
ameter D25 and is centered on the R band center of 
light of the galaxy. Bottom left panel: The synthe- 
sized H I line profile obtained by integrating the data 
cube along both spatial axes. Bottom right panel: 
The rotation Vrot sin i and dynamical major axis po- 
sition angle F as a function of radius in arcseconds, 
obtained from the tilted ring model. The data for the 
approaching side (circles; solid lines) and the reced- 
ing side (triangles; dashed lines) are presented super- 
posed. 



Fig. 4.— VLA H I synthesis results for UGC 3685. 

Top left panel: H I cohimn density contours detected 
at the VLA, with light gray contour at 1.1 x 10^^ cm~^ 
and solid contours at {3.3,13.2,23.1,33.0,...} x 
lO^'' cm^^, overlaid on the R band image obtained by 
KHL. The first solid contour represents the 3a noise 
level. Top right panel: The same H I column den- 
sity contours as in the top left panel, plotted over a 
linear grayscale column density map. The circle is of 
diameter D25 and is centered on the i?-band center of 
light of the galaxy. Center left panel: The position- 
velocity diagram obtained from the H I data cube. 
The slice was made along the dynamical major axis 
F = 299° determined by averaging the dynamical po- 
sition angle over radii r < i?25- The x-axis indicates 
arcminutes northwest of the dynamical center, and 
the dotted vertical lines identify i?25- Center right 
panel: The H I observed velocity field with contours 
and grayscales from 1730 km s~^ (east, light shading) 
to 1830 km s""'^ (west, dark shading) in increments of 
10 km s-^ The 1780, 1810, and 1830 km s'^ con- 
tours are identified in the panel. The circle is of di- 
ameter D25 and is centered on the J?- band center of 
light of the galaxy. Bottom left panel: The synthe- 
sized H I line profile obtained by integrating the data 
cube along both spatial axes. Bottom right panel: 
The rotation Vrot sin i and dynamical major axis po- 
sition angle F as a function of radius in arcseconds, 
obtained from the tilted ring model. The data for the 
approaching side (circles; solid lines) and the reced- 
ing side (triangles; dashed lines) are presented super- 
posed. 
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Fig. 5. — VLA H I synthesis results for 
NGC 3596. Top left panel: H I column density 
contours detected at the VLA, with contours at 
{4.05,16.2,28.35,40.5,...} x lO^^ cm-^, overlaid on 
the i?-band image obtained by KHL. The first con- 
tour represents the 3(t noise level. Top right panel: 
The same H I column density contours as in the top 
left panel, plotted over a linear grayscale column den- 
sity map. The circle is of diameter D25 and is centered 
on the i?-band center of light of the galaxy. Cen- 
ter left panel: The position- velocity diagram obtained 
from the H I data cube. The slice was made along the 
dynamical major axis T = 80° determined by averag- 
ing the dynamical position angle over radii r < i?25- 
The a;-axis indicates arcminutcs west of the dynami- 
cal center, and the dotted vertical lines identify i?25- 
Center right panel: The H I observed velocity field 
with contours and grayscales from 1130 km s~^ (west, 
light shading) to 1260 km s~^ (east, dark shading) 
in increments of 10 km s'^ The 1170, 1200, and 
1220 km s~^ contours arc identified in the panel. The 
circle is of diameter D25 and is centered on the R- 
band center of light of the galaxy. Bottom left panel: 
The synthesized H I line profile obtained by integrat- 
ing the data cube along both spatial axes. Bottom 
right panel: The rotation Vrotsini and dynamical 
major axis position angle F as a function of radius 
in arcseconds, obtained from the tilted-ring model. 
The data for the approaching side (circles; solid lines) 
and the receding side (triangles; dashed lines) are pre- 
sented superposed. 



Fig. 6. VLA H I synthesis rcsuhs for UGC 6429. 
Top left panel: H I column density contours detected 
at the VLA, with contours at {5.4, 16.2, 27, 37.8, . . .} x 
10"'^^ cm~^, overlaid on the R band image obtained by 
KHL. The first contour represents the 3a noise level. 
Top right panel: The same H I column density con- 
tours as in the top left panel, plotted over a linear 
grayscale column density map. The circle is of diam- 
eter D25 and is centered on the i?-band center of light 
of the galaxy. Center left panel: The position- velocity 
diagram obtained from the H I data cube. The slice 
was made along the dynamical major axis F = — 9° 
determined by averaging the dynamical position angle 
over radii r < i?25- The x-axis indicates arcminutes 
north of the dynamical center, and the dotted vertical 
lines identify R25. Center right panel: The H I ob- 
served velocity field with contours and grayscales from 
3700 km (south, light shading) to 3770 km s'^ 
(north, dark shading) in increments of 5 km s"-'^. The 
3710, 3735, and 3750 km s~^ contours are identified 
in the panel. The circle is of diameter £'25 and is 
centered on the R band center of light of the galaxy. 
Bottom left panel: The synthesized H I line profile ob- 
tained by integrating the data cube along both spatial 
axes. Bottom right panel: The rotation Vrot sin i and 
dynamical major axis position angle F as a function 
of radius in arcseconds, obtained from the tilted-ring 
model. The data for the approaching side (circles; 
solid lines) and the receding side (triangles; dashed 
lines) are presented superposed. 
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Fig. 7. — VLA H I synthesis results for 
NGC 4688. Top left panel: H I column density 
contours detected at the VLA, with contours at 
{7.2, 21.6, 36.0, 50.4, . . .} x 10^^ cm~^ overlaid on the 
i?-band image obtained by KHL. The first contour 
represents the 3cr noise level. Top right panel: The 
same H I column density contours as in the top left 
panel, plotted over a linear grayscale column density 
map. The circle is of diameter D25 and is centered on 
the iZ-band center of light of the galaxy. Center left 
panel: The position-velocity diagram obtained from 
the H I data cube. The slice was made along the 
dynamical major axis F = 3° determined by averag- 
ing the dynamical position angle over radii r < i?25- 
The X-axis indicates arcminutes north of the dynam- 
ical center, and the dotted vertical lines identify i?25- 
Center right panel: The H I observed velocity field 
with contours and grayscales from 960 km s^^ (south, 
light shading) to 1010 km s^^ (north, dark shad- 
ing) in increments of 5 km s^^. The 960, 980, and 
1000 km s^^ contours are identified in the panel. The 
circle is of diameter D25 and is centered on the R~ 
band center of light of the galaxy. Bottom left panel: 
The synthesized H I line profile obtained by integrat- 
ing the data cube along both spatial axes in the re- 
gion containing flux from NGC 4688. The profile does 
not include flux from the companion. Bottom right 
panel: The rotation Vrotsini and dynamical major 
axis position angle F as a function of radius in arcsec- 
onds, obtained from the tilted-ring model. The data 
for the approaching side (circles; solid lines) and the 
receding side (triangles; dashed lines) are presented 
superposed. 



Fig. 8.— Wide field VLA H I synthesis data for 
CGCG 043-029, the companion of NGC 4688. Top 
left panel: The H I column density of NGC 4688 and 
companion. The contours and grayscales are as in fig- 
ure 1^. Top right panel: The H I column density of the 
companion object alone, overlaid on the POSS image 
of the region. Bottom left panel: The integrated line 
profile of the companion alone, obtained by integrat- 
ing both spatial axes of the data cube in the region 
containing flux from the companion. Bottom right 
panel: the POSS image of CGCG 043-029 without 
the H I contours, showing clearly the two distinct op- 
tical maxima. 



Fig. 9.— VLA H I synthesis resuhs for NGC 5474 
obtained from VLA archive data originally reported 
by RDH. Top left panel: H I column density 
contours detected at the VLA, with contours at 
{1.8, 10.8, 19.8, 28.8, . . .} x 10^^ cm^^ overlaid on the 
iZ-band image obtained by KHL. The first contour 
represents the 3cr noise level. Top right panel: The 
same H I column density contours as in the top left 
panel, plotted over a linear grayscale column density 
map. The circle is of diameter D25 and is centered 
on the JZ-band center of light of the galaxy. Center 
left panel: The synthesized H I line profile obtained 
by integrating the data cube along both spatial axes. 
Center right panel: The H I observed velocity field 
with contours and grayscales from 259 km s~^ (west, 
light shading) to 297 km s~^ (east, dark shading) in 
increments of 4 km s~^. The 261, and 295 km s"^ 
contours are identified in the panel. The circle is 
of diameter D25 and is centered on the dynamical 
center of the galaxy as determined by tilted-ring fit- 
ting. Bottom left panel: The rotation Vrot sin i and 
dynamical major axis position angle F as a function 
of radius in arcseconds, obtained from the tilted-ring 
model. The data for the approaching side (circles; 
sofid fines) and the receding side (triangles; dashed 
lines) are presented superposed. 
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Fig. 10.— VLA H I synthesis results for NGC 5701 
obtained from VLA archive data originally observed 
by S. E. Schneider. Top left panel: H I column den- 
sity contours detected at the VLA, with contours at 
{2.9, 7.8, 12.7, 17.6, . . .} x lO^^ cm'^, overlaid on the 
i?-band image obtained by KHL. The first contour 
represents the 3ct noise level. Top right panel: The 
same H I column density contours as in the top left 
panel, plotted over a linear grayscale column density 
map. The circle is of diameter D25 and is centered on 
the i?-band center of light of the galaxy. Center left 
panel: The synthesized H I line profile obtained by in- 
tegrating the data cube along both spatial axes. Cen- 
ter right panel: The H I observed velocity field with 
contours and grayscales from 1440 km s"""^ (west, light 
shading) to 1560 km s~^ (east, dark shading) in incre- 
ments of 10 km s-^ The 1450, 1470, and 1550 km s'^ 
contours are identified in the panel. The circle is of 
diameter 1)25 and is centered on the i?-band center of 
light of the galaxy. Bottom left panel: The rotation 
Vrot sin? and dynamical major axis position angle F as 
a function of radius in arcseconds, obtained from the 
tilted-ring model. The data for the approaching side 
(circles; solid lines) and the receding side (triangles; 
dashed lines) are presented superposed. 



Fig. 11.— VLA H I synthesis resuUs for UGC 12732 
obtained from VLA archive data originally observed 
by Schulman et al. (1997). Top left panel: H I col- 
umn density contours detected at the VLA, with light 
gray contour at 6.7 x 10^* cm~^ and solid contours 
at {2.0, 11.4, 20.8, 30.2, . . .} x lO^^ cm'^, overlaid on 
the i?-band image obtained by KHL. The first solid 
contour represents the 3a noise level. Top right panel: 
The same H I column density contours as in the top 
left panel, plotted over a linear grayscale column den- 
sity map. The circle is of diameter D25 and is centered 
on the i?-band center of light of the galaxy. Center 
left panel: The synthesized H I line profile obtained 
by integrating the data cube along both spatial axes. 
Center right panel: The H I observed velocity field 
with contours and grayscales from 680 km s~^ (south, 
light shading) to 810 km s~^ (east, dark shading) in 
increments of 5 km s~^. The 720, 760, and 810 km s~^ 
contours are identified in the panel. The circle is of 
diameter D25 and is centered on the i?-band center of 
light of the galaxy. Bottom left panel: The rotation 
Vrot sin« and dynamical major axis position angle F as 
a function of radius in arcseconds, obtained from the 
tilted-ring model. The data for the approaching side 
(circles; soHd lines) and the receding side (triangles; 
dashed lines) are presented superposed. 

Fig. 12. — Illustrations of the observed relation- 
ships between the morphological asymmetry parame- 
ters (Af^Hi and ^/,fl) and the dynamical parameters 
measuring position angle change with radius AF(r), 
position angle change from approaching to receding 
side AT{9), K-otsini(r) curve asymmetry S2, and in- 
tegrated H I line profile asymmetry A„ . 

Fig. 13. — Illustrations of the relationships between 
H I morphology and the global H I parameters, from 
top to bottom, mass-to-light ratio in solar units, H I 
mass, H I diameter in terms of optical diameter -D25, 
and projected distance to the nearest neighbor. The 
positions of UGC 3685 and NGC 3596, which have no 
neighbors within l°are indicated by vertical arrows. 



Fig. 14. — Illustrations of the relationships between 
rotation curve asymmetry S2 and the global H I pa- 
rameters examined in Figure 
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